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Phase-sensitive properties of light play a crucial role in a variety of quantum optical phenomena, 
which have been mostly discussed in the framework of photoelectric detection theory. However, 
modern detection schemes, such as arrays of on-off detectors, are not based on photoelectric counting. 

We demonstrate that the theory of homodyning with such click counting detectors can be established 
by using a proper detection model. For practical applications, a variety of typically occurring 
imperfections are rigorously analyzed and directly observable nonclassicality criteria are studied. 
Fundamental examples demonstrate the general functionality of our technique. Thus, our approach 
of homodyne detection with on-off detector systems is able to bridge the gap between imperfect 
detection and the phase resolution demands for modern applications of quantum light. 

PACS numbers: 42.50.Ar, 03.65.Wj, 85.60.Gz 


I. INTRODUCTION 

Phase-sensitive measurements have a fundamental im¬ 
pact on uncovering quantum properties of radiation 
fields [IHS] . Applications of quantum effects can be found 
in the vast fields of quantum information science and 
quantum metrology [IHZl- Those are typically studied in 
terms of homodyne measurements with high intensities 
and the photoelectric detection theory BM- However, 
applications of light quanta in the single- or few-photon 
domain faces some flaws as realistic photon-number re¬ 
solving detectors are often not available pT)] . 

The nonclassicality of light is usually determined with 
reference to classical coherent states mnn]. More pre¬ 
cisely, any quantum state can be represented in the ba¬ 
sis of these classical states by the Glauber-Sudarshan 
P phase-space representation Whenever the P 

function is nonnegative, it can be interpreted as a prob¬ 
ability distribution of a classical ensemble of electromag¬ 
netic waves. If the P function exhibits negativities, no 
such interpretation can be made, and the state of light 
is therefore a nonclassical one. Hence, the categorization 
in classical or nonclassical radiation fields plays a crucial 
role. 

In particular, phase dependent quantum phenomena 
are often studied using interferometric measurement 
schemes such as homodyne detection [IMH]- In a homo¬ 
dyning setup, a signal (SI) is superimposed on a beam 
splitter with a local oscillator (TO), with a controllable 
phase. This is usually done with a strong LO and formu¬ 
lated in terms of the photoelectric detection theory. Rig¬ 
orous analysis for low LO intensities turns out to be more 
involved but can be analyzed and implemented [I^21j . 
In addition, homodyne detection serves as a foundation 
for quantum state reconstruction techniques [mill] for 
uncovering all nonclassical features of light. 

Hand in hand with the quantum nature of the photon 
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comes its particle nature, which may be confirmed by an- 
tibunching effects [2^25] . Thus, to gain further knowl¬ 
edge of the nature of light, the generation and analysis of 
states in a few-photon regime has gained more and more 
importance This has brought up the necessity for 

adequate detectors. For example, hybrid detectors [28] . 
superconducting detectors [11130], and click-counting de¬ 
tectors [3TH50] , consisting of multiple avalanche photodi¬ 
odes (APDs) in the Geiger mode, may meet this demand. 
The latter ones offer a nonlinear, but well-defined, statis¬ 
tics [371139] . Quantum correlations in the few photon 
domain can be inferred with such click detectors [351110]- 
m], even without the need for additional data process¬ 
ing [45]. To get the full quantum picture, it is desirable to 
investigate both particle- and phase-dependent phenom¬ 
ena simultaneously. In Ref. [10] , a step towards a theory 
of phase-resolving click counting was made. We aim at 
extending this investigation by giving a more general ap¬ 
proach and lay the foundation for further applications, 
such as quantum state reconstruction 113 - 

In the current contribution, we formulate the theory 
for general homodyne detection measurements with click 
detectors. This includes balanced and unbalanced de¬ 
tection with four or more port homodyning setups. The 
verification of nonclassical features is formulated in terms 
of measured second order correlations, and the straight¬ 
forward generalization to higher-order correlations is out¬ 
lined. We rigorously investigate the influence of imperfec¬ 
tions in such measurement scenarios, which includes, e.g., 
the impact of detector imperfections as well as LO fluc¬ 
tuations and mode mismatch between LO and SI. More¬ 
over, the identification of quantum correlations between 
multiple SI fields is shown using several homodyne detec¬ 
tor settings. 

This paper is structured as follows. In Sec. [H] we give 
an overview of the treatment of click counting detec¬ 
tors, reviewing results achieved on this topic so far. The 
click detection in general four port homodyning schemes 
is derived in Sec. |HI[ including the cases of unbalanced 
and balanced scenarios. The influence of imperfections 
is studied in Sec. m Multi-port homodyne detection 
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schemes are treated in Sec. [V| by means of the example 
of a balanced eight-port setup. The discourse on phase- 
sensitive click counting is then completed in Sec. lYll 
where we discuss homodyne correlation measurements 
between multiple signal fields. Finally, we summarize 
and conclude in Sec. m 


II. CLICK-COUNTING THEORY 


A. Click counting versus photoelectric counting 


Most commonly, when detection processes of light are 
discussed in quantum optics, the theory of photoelectric 
measurement is examined. For the single mode case, the 
resulting photoelectric counting statistics reads m 


Pn = 


n! 


( 1 ) 


where : • • •: is the normal ordering prescription, 77 is the 
quantum efficiency, and v is the dark count rate |48j . 
The photoelectric statistics is (i) the true photon number 
statistics for a perfect detection scenario, ry = 1 and iz = 
0, and (ii) a true Poisson statistics for coherent light. 




CC 
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FIG. 1: (Color online) The top illustration shows the in¬ 
tegrated symbol “CC”, which we will further use for click 
counting detectors of any type. Two realizations of click de¬ 
tectors with A = 4 APDs are additionally shown. In the 
spatial multiplexing setup (middle scheme) the incident light 
is equally divided by multiple 50:50 beam splitters. Time-bin 
multiplexing is an equivalent realization. In the detector ar¬ 
ray scenario (bottom scheme), an array of APDs is equally 
illuminated. 

Alternatively, detector systems based on APDs can be 
applied. In such a configuration, each APD acts as an 
on-off detector Ea mi- Realizations of detector sys¬ 
tems consisting of multiple APDs are detector arrays 
or multiplexing setups (cf. Fig. [^. Time-bin multi¬ 
plexing is related to the spatial multiplexing setup (cf., 
e.g., EH EZl mi US])- Experimental characterizations 


of such click detector systems have been performed in 

Refs. [501 ED- 

A fundamental approach to retrieve a kind of photon 
number resolution is to equally distribute the incident 
light onto N APDs. This yields a total number of k 
clicks, 0 < k < N. The measured statistics of such click 
detector devices can be described with a quantum version 
of a binomial statistics |M] : 

where N is the number of APDs and Ck gives the proba¬ 
bility to measure k clicks. The mean value of this statis¬ 
tics may be given by the expectation value of the operator 
tt: 

N 

(tt) = ^ fccfc, with TT = A . (3) 

k=0 

This operator is a nonlinear function of the photon num¬ 
ber operator. Note that we use a rescaled version (factor 
N) of the operator introduced in Ref. [52] . 

Due to the finite nature of the click-counting statistics, 
0 < /c < A, all possible events are accessible in experi¬ 
ments. In contrast, the photoelectric detection model Q 
(only in theory) allows one to have arbitrarily high n val¬ 
ues. For instance, the prediction of a coherent state |q;) 
yields Pn ^ 0 for all n values, which cannot be achieved 
with a finite number of measured events. 

Let us give another example to underline the differ¬ 
ence between click counting and photoelectric counting 
theory. Three prominent states will serve as our test 
states throughout this paper. One is a classical coher¬ 
ent state |a), and, to illustrate the basic differences from 
nonclassical light fields, we analyze the even (-I-) and odd 
(—) coherent states [55] , 

l+a)±|-a) 

\(^±/ = y ( 4 ) 

y/ 2 [I±exp(- 2 |a| 2 )] 

In Fig. we compare the click statistics [(al)-(a3)] 
and photoelectric counting statistics [(bl)-(b3)] of these 
states for a perfect detection scenario, rj = 1 and 1 / = 0. 
Let us stress that the click statistics is limited to the 
displayed plot range, 0 < k < 8 = N, whereas the 
photoelectric statistics is not. Figures [^al) and ibl), 
|^a 2 ) and[^b 2 ), and [^a3) and ^b3) show the coher¬ 
ent, even coherent, and odd coherent states, respectively. 
It is clearly visible that the true photon statistics of the 
(even)odd coherent state includes only (even)odd photon 
numbers. This feature is not distinct in the correspond¬ 
ing click statistics, even for the studied ideal detection. 


B. Moment based nonclassicality probes 

Based on the variance of the click statistics, it has 
been demonstrated that it is possible to infer nonclas¬ 
sical features of quantum light, yielding the notion of 
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second-order moment based constraint for classical light 
is 


0<(:[Ad]2:) =detM(2) (8) 

= ((:!:) oJ'o) = (J (:(Ad)2:)) 

= '^''^((:(7V-7r):) (KiV-d)^) 

If, in addition, we consider two click-counting detectors, 
we can also formulate cross-correlation criteria by the 
corresponding minors [52] . For example, the second order 
criterion is 

0 < (:(Adi)2:)(:(Ad2)2:) - (iAttiAT r^:)^, (9) 

where each click detector {i = 1,2) is characterized by the 
operator tti [cf. Eq. (|^]. These cross correlations have 
been recently used to experimentally uncover quantum 
correlations between two light beams |45j . 

Finally, it is worth mentioning that the moments 
for any number of click-counting devices can be di¬ 
rectly obtained from the measured joint click statistics 
Cfci,fc 2 ,... |52) . This yields joint moments of the form 


FIG. 2: (Color online) (al)-(a3) The click statistics as given 
in Eq. ® for A = 8 is compared to (bl)-(b3) the corre¬ 
sponding photoelectric counting statistics in Eq. 0- The 
considered states are (top) the coherent state |a), as well as 
(middle) odd and (bottom) even coherent states For all 
plots we take a = 2, rj — 1, and = 0. 

sub-binomial light [531 Hi] • In general, it was also shown 
that higher-order ordered moment criteria can be used 
to asses these criteria |53|. For this reason the matrix of 
Ath-order moments of the click statistics was defined as 

, (5) 

V / m,m '—0 

where the even integer K satisfies K < N. If we do not 
restrict ourselves a particular moment order, we skip the 
upper index and write M . For classical radiation fields 
this matrix is non-negative, M > 0, and a violation of 
this property uncovers quantum light. 

It might be more convenient to formulate the modified 
matrix of moments in terms of 

{■.{N - 7r)'"+™ :) = (: :) (6) 

or central moments 

(:(Ad)™+™';) = (:(7r - (:7r:))'"+™':). (7) 


Ni,N 2 ,... (ki\ 

tn*r^)=n"r i: 

i i ki—m\,k2—m2,... i \miJ 

( 10 ) 

which makes the matrix of click-counting moments ap¬ 
proach easily accessible. Based on these moments, a con¬ 
struction of higher-order or cross-correlation nonclassi- 
cality conditions can be formulated for any number of 
click detector systems; see also (52). 

C. Outline 

In summary, the click counting approach is a directly 
applicable technique to infer phase-insensitive quantum 
properties of radiation. The click statistics (pj) itself can 
be obtained from the photoelectric statistics TB only via 
a non bijective mapping from the infinite number of pho¬ 
ton numbers to the finite spectrum of click counts (cf. 
Ref. 1331) • Due to this feature, it is evident that the im¬ 
plementation of click detectors, in the theory of phase- 
sensitive measurement, has to be carried out with great 
care. In the same way, nonclassicality probes have to 
be adjusted to avoid misleading interpretations of the 
probed state of light. In the remainder of this work, we 
will perform such an analysis for a number of established 
homodyning schemes in quantum optics. 


The resulting ways for determining nonclassicality are 
identical for all the different formulations of matrices of 
click counting moments, which is proven in Appendix [A] 
The non-negativity of all these matrices of moments can 
be inferred from their principal minors. For example, a 


III. FOUR PORT HOMODYNE DETECTION 

We will start our treatment of phase sensitive measure¬ 
ment using click detectors with the most prominent setup 
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in homodyne detection: the four port scheme (see Fig.|^. 
Applications of this measurement scheme, in terms of 
photoelectric detection theory, include fundamental ap¬ 
plications such as quantum state tomography liiniiia. 
In this measurement scheme, a SI is superimposed with 
the LO on a beam splitter. The latter is described by 
a transmission coefficient t and a reflection coefficient r. 
Let us emphasize that we will describe this setup em¬ 
ploying click counting detectors instead of the standard 
scenario with photoelectric detectors. 


r^cc 



FIG. 3: (Color online) The homodyne four-port scheme. The 
LO, having a controllable phase, and the SI are superimposed 
on a |tp:|r|^ beam splitter. The output beams 1 and 2 are 
each detected with a click detector. 


The input-output relation of the beam splitter can be 
expressed by a unitary transformation matrix, 



where -I- |rp = 1. The click detectors in positions 1 
and 2 in Fig. yield a joint click-counting statistics 



( 12 ) 


where hi = a]hi is the photon number operator, rji is 
the quantum efficiency, t'i is the dark count rate, and Ni 
is the number of APDs of the click detector at position 
i = 1,2. The reference beam is a coherent light field, 
\P) = |/d)LO- Thus, we can insert the photon number 
operators in Eq. (12 1 in the form 



Note that this form corresponds, up to a scaling, to dis¬ 
placed photon number operators. 


^(t) = (asi - 7)^(asi - 7 ): (15) 


where h = dgjdsi is the photon number operator of the 
SI. Eventually, the full click-counting statistics of a four 
port detection scheme is 


Cfei,fc2 — 


^g-r;i|t|^ri(-r/3/t)/Afi-!/i^ 


ATi-fci 


N 2 

k2 


)(' 


-'^2k|^n(i*/3/r*)/Ar2-i'2 


N2-k2 




This joint click statistics of the general four port ho¬ 
modyne detector may be used to infer nonclassical light 
fields via moment criteria. For example, the normally or¬ 
dered click-counting variances or cross correlations certify 
quantumness if 

0 >(:[A7ri]2:) = A^2g-2i2i (■^7) 

0 >(:[A7f2]2:) = (ig) 

0 >(:[A7ri]2:)(:[A7r2]2:) - (iAttiA 7r2:)^ (19) 


In the following, we will focus on two specific four port 
homodyning schemes with some relevance in quantum 
optics. 


A. Unbalanced detection 


Let us first consider the unbalanced measurement |56] 
with the click detector only in channel 1 of Fig. That is, 
contributions of the detector at position 2 are traced out. 
Alternatively, this corresponds to the case 1^2 = V 2 = 0, 
with a joint click counting statistics ( 121 , which has the 


property that = 0 for k 2 yf 0. For convenience, we 
may replace the notations Ck,o = Ck, Ni = N, r]i = rj, 
and i 2 i = v in this case. Thus, for a LO |/3), we get 




N-k 




using the beam splitter transformation ( 11 ), with 774 = 
|tp 77 being an overall quantum efficiency and 7 = —r/3/t. 
Recently, it has been shown in Ref. [47] that such unbal¬ 
anced homodyne setups reveal nonclassical features in 
terms of click counterparts of s-parametrized quasiprob¬ 
abilities. 

If we decompose the coherent displacement in the form 
7 = 1716 *“^, the phase-sensitive mean click counts are 
given as the expectation value of the operator 'iT((p) [cf. 
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Eq. (|^]. For instance, we get moments of the form 

{:[N - 7r(vj)]'”:) = (21) 

Now, nonclassicality criteria can be constructed. For ex¬ 
ample, we have nonclassical light for the LO phase (p if 


0 > (: [ATf{ip)Y :) = N^e 


2„-2u / 


Ae-t"(7) 


1 2 


:). ( 22 ) 


Higher-order nonclassicality can be identified with 
the matrix of phase-sensitive click-counting moments, 
= ((:7r'"+'” {p)-))m,m', which is non-negative for 
classical light, M{ip) > 0. For example, the fourth order 
nonclassicality condition is 


( 1 

0 >det ^((^)=det I (: 7 r((y 9 ):) (: 7 r^((/ 3 ):) 

\(:7r2(^):) (:,r3(^):) 


{:nHpy.)) 

(23) 


Examples of these phase-sensitive variances are given 
in Fig. for the even and odd coherent states in Eq. Q. 
A phase-sensitive verification of nonclassicality can be 
observed for the even coherent state, while the nonclas¬ 
sicality of the odd coherent state is not revealed by this 
particular nonclassicality probe. The effective quantum 
efficiency in this unbalanced homodyne detection setup 
is chosen to be T]t = 32%. For the plot of quantity (22 1 , 
the analytic result of the expectation value, 


(Q!±|:exp [-A(a - 7 )’l'(a - 7 )] :|a±) 

-AKI^ cosh[ 2 ARe( 7 *a)] 

® - l±e-2|«P 

cos[ 2 AIm( 7 *a)] 


(24) 


± e 


1 ± e-2|“l" 


was employed. It also yields the analytical results for the 
following discussions. 

Let us additionally study the verification of nonclassi¬ 
cality for different N values. Click detection schemes in 
Fig.[T] may also employ CCD cameras with a high num¬ 
ber of APDs; see, e.g., [STllSH] for theoretical and exper¬ 
imental studies using CCD cameras. In Fig. we show 
the N dependence of the verified nonclassicality for an 
even coherent state for second- and fourth-order criteria 
[see Eqs. (|^ and (23), respectively]. The negativities 


are typically more pronounced for the fourth order cri¬ 
terion (right plot) in comparison with the second order 
one (left plot). It can be observed that higher numbers 
of APDs are advantageous for those criteria too. In the 
limit N ^ 00 , the negativities converge to those val¬ 
ues which are expected for an unbalanced homodyning 
scheme with a photoelectric detector. 


B. Balanced detection 

Now we will study the balanced case of the four-port 
scheme in Fig. That is, the beam-splitter transforma¬ 
tion in Eq. ([IT|) is specified hy t = r = 1/ \f2. In addition. 



0 ;r 2;r 


FIG. 4: (Color online) The application of nonclassicality 
criterion in Eq. \22\ is shown as a function of the phase p 
for !«-) (dashed line) and |q;+) (solid line), where a = 1 
and \j3\ = 4. Negativity of the normally ordered variance 
determines nonclassicality. The click detection parameters 
are 77 = 50%, 72 = 0, A = 8, and the beam splitter in Eq. (11) 
is characterized by t = 4/5 and r = 3/5. 



FIG. 5: (Color online) (left) The application of the second- 
order nonclassicality criterion in Eq. ( |22[ ) for a fixed phase 
(p = 7 r /2 and (right) the fourth-order criterion in Eq. (23) 
are shown as a function of the number N of APDs for the 
even coherent state |a+) with a = 1. The displayed values 
determine nonclassicality, and the insets illustrate the behav¬ 
ior for small N values. The click detection parameters are 
7 = 50%, V = Q, |/3| = 4, and the beam splitter in Eq. (Ill is 
characterized hy t = 4/5 and r = 3/5. 


we assume that both detectors have the same character¬ 
istics, i.e., Ni = N 2 = N, rji = ri 2 = 7 , and vi = V 2 = v. 
On this basis, the joint click-counting statistics in (12) 
reads 



(25) 


The theory of balanced homodyne detection with click 
detectors has been established in |46] . Here, we will recall 
some of its features and complete the discussion of some 
aspects not considered previously. 

When working with photoelectric detectors in the bal¬ 
anced four-port scheme, one subtracts the photoelectric 
counts of the two detectors from one another. In the limit 
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of strong LO, one measures the quadrature 

x{(p) = + dgie*”^ (26) 

of the signal. Analogously, the click difference counts can 
be analyzed, which yields the moments of a nonlinear 
click-quadrature operator [46], 

(:A™((p):)=(:(7ri-7r2)-:) (27) 

j=0 ^ ^ 

with 

X{ip) = 2Afe“2Nl/5p-'";g-2w”-sinh 

The quadrature x{ip) appears here in a hyperbolic sine 
function, which is a skew-symmetric, analytical function, 
sinh(a;) = {e^ — e“^)/2 = — sinh(—a:). It transforms the 
features we expect from a quadrature and thereby justi¬ 
fies the notion of a nonlinear quadrature. Consequently, 
a nonlinear squeezing condition can be given as 

0>(:[AA(<^)]2:). (29) 

Moreover, additional phase-sensitive information can 
be obtained from the jointly measured click statistics, as 
it reveals some fundamental differences from the photo¬ 
electric detection model. In particular, the sum of click 
counts does, in contrast to its photoelectric counterpart, 
show phase-dependent behavior. The number of the sum 
of photoelectric counts is independent of the phase, 

{rjhi -\-v) + (77712 + v) = Tydgjdsi + 'q\f}\'^ + 2v. (30) 

The sum of click counts, 

7fi-|-7r2 = 2A—“'^:e“^"cosh 

(31) 

clearly shows a phase-dependent behavior. Thus, a non- 
classicality criterion based on the variance of the sum can 
be formulated similarly to the case of the click difference. 

Examples of both criteria are shown for the even 
and odd coherent states 0 in Fig. Nonlinear 

squeezing can be observed for the even coherent state, 
whereas a sub-shot-noise variance of the sum of clicks, 
(:[A(7ri-|-^2)]^:) < 0, is visible for the odd state. More¬ 
over, it is worth mentioning that we are not limited to 
weak or strong LOs. 


^1^1 ^ 


rm, 

2N‘ 




:. (28) 


IV. IMPERFECTIONS 



FIG. 6 : (Color online) The result of phase-dependent 

click statistics for a balanced homodyning scheme. The left 
plot shows the application of nonlinear squeezing criterion in 
Eq. ( |29| l, and the right plot (scaled xlO^) depicts the corre¬ 
sponding variance of the sum of clicks. The plotted states 
are |a-) (dashed line) and |a+) (solid line), where q = 1 
and \j3\ = 4. The detection parameters are 77 = 50%, 1 ^ = 0, 
N = 8. 


click operators, e.g., Eqs. 0 and it is directly evi¬ 
dent that click detectors respond differently to attenua¬ 
tions than photoelectric detectors. Moreover, the impact 
of a mode mismatch between the LO and the SI has not 
been studied so far. We will therefore study these realis¬ 
tic perturbations for the determination of nonclassicality 
while restricting ourselves to the fundamental case of un¬ 
balanced homodyne detection. 


A. Imperfections of the click-counting detectors 


One imperfection is due to the dark count rate iz [see 
Eqs. 0 and ([^]. That is, some clicks are recorded 
even if there was no SI or LO field. Let us consider the 
nonclassicality criteria, i.e., moments of the form (211. 
Comparing the case zz = 0 and a nonzero zz, we can de¬ 
compose the matrix of moments as 


/ / \ 1^/2J 

\ / m,m —0 

/ r /I TJ+ - / X / \ L-^/2j 

= (e-^['"+'"](:[Ae-^”(T')]'"+'" :)) (32) 

V / m,m'=0 

=T, with n = diag (e-°",..., e'Lf . 


This transformation property between the cases with and 
without dark counts allows us to state the following: If 
nonclassicality can be detected for no dark counts, zz = 0 , 
then nonclassicality can be detected for a finite dark 
count rate, zz > 0. Thus, the impact of dark counts, 
excluding the case zz —> 00 , is not an issue for the verifi¬ 
cation of nonclassicality. It solely scales the actual values 
of the minors, e.g., 

X (33) 


Let us now discuss the influence of imperfections. Thus, we can assume for our theoretical studies zz = 0. 
which is crucial for the verification of nonclassical- In experiments zz can also be estimated, e.g., zz « 0 [50] 
ity [451 [551 - I5T] . From the nonlinear structure of the mean or zz « 0.5 [15] . 
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Such a simple treatment is not possible when consid¬ 
ering a non unit quantum efficiency, 77 < 1. This analysis 
has to be performed specifically for the desired target 
state. Such a test, in advance of an experiment, is help¬ 
ful for estimating the overall efficiency one requires to 
infer nonclassicality. 

Let us take a closer look at the even coherent state (see 
Fig.0. For this state the appearance of nonclassicality 
seems to be quite independent of the quantum efficiency. 
That is, the value of the negativity is diminished for low¬ 
ered efficiency, but the negativity is present for the 7r/2 
and 37 r /2 phases for any 77 > 0. However, the signifi¬ 
cance, with which the nonclassicality might be verified 
in experiments, decreases. 


-<:[A7r(y)]2:> 



FIG. 7: (Color online) The nonclassicality criterion \22\ is 
shown depending on the quantum efficiency 77 and the phase 
ip. For better visibility, the negative variance, —(:[A'n'(i^)]^:), 
is plotted. We consider an even coherent state |a+}, where 
0 = 1 and /3 = 4, the click detector parameters are v = Q and 
N = 8, and the beam splitter is characterized by t = 4/5 and 
r = 3/5. 


B. Saturation effects 

The number N of APDs is finite. It can be shown 
(cf. [311 Sn]) that we approach a true displaced number 
operator, 

= (:iV™(l-e-^"(^))-:) « (:[77*n(7)]'":), 

(34) 

in the limit iV —> 00 . Another version of the limit would 
be that the powers of photon numbers are comparably 
small, (:n'"( 7 ):) <C IV™, for ttt- S N \ {0}. 

Both approximations are typically not justified. Even 
worse, accepting these approximations would yield non¬ 
classicality for classical coherent states [M]. To counter 
these fake effects, all our nonclassicality criteria are for¬ 
mulated in terms of moments of the click statistics. This 
allows us to treat all ranges of intensities without ap¬ 
proximations. This also means that saturation effects are 


included in our click-counting theory. Namely, the case 
(:h'"( 7 ):) > IV™ is properly described. For tti = 1 this 
means that the mean photon number ( 77 ,( 7 )) can be on 
the same order as or exceed the number of on-off diodes, 

N. 

Let us show that fake nonclassicality will not occur 
for high intensities. Thus, we may assume a quantum 
SI with many photons and a classical SI with the same 
intensity. Basically, the high-intensity limit of both cases 
yields that all APDs click at the same time, i.e.. 


Jo for fc = 0,..., — 1, 

1 for k = N. 


(35) 


Note that the same is true if the LO is strong. For the 
moments [see formula ( 10 )] we get 


(:7r™((p):) = IV" 


N 

E 


^ r) 

—m \m/ 


N” 


(36) 


which consistently yields a positive matrix of moments, 
(^]\[m+m > 0 , for the classical and nonclassical 

states. In conclusion, (i) too high intensities are not 
helpful for the determination of nonclassicality, and (ii) 
a proper detector description will not yield fake nonclas¬ 
sicality, even in case the of saturation. 


C. Imperfections due to the local oscillator 

Imperfections in homodyne detection may also stem 
solely from the LO. Typical examples in realistic mea¬ 
surement scenarios are fluctuations of the LO itself or a 
mode mismatch between the LO and the SI fields. Let us 
discuss the implications for such sources of errors. Explic¬ 
itly, the impact on 777 th order moments will be studied, 

(:[Af-7r((/7)]™:) = (iVe-'^)™(:e-T^''‘(“si-7)t0si-7).)^ 

(37) 

with 7 = —r/3/t (ip = arg 7 ) being a linear function of 
the amplitude of the LO state, |/3)lo- 

For small amplitudes some noise of the LO may occur, 
e.g., due to thermal fluctuations. For the balanced four- 
port homodyne detection, phase and amplitude noise 
have been studied in Ref. |46]. Here, we will focus on 
perturbations in the unbalanced scenario. If the source 
of noise is a classical one, we may describe this effect by 
a classical probability distribution, Phoil) > 0 , of the 
(scaled) LO amplitude. Thus, we get a convolution of 
the moment with noise: 

(:[IV-7r(:^)]™:)p,, (38) 

= J dV-PLo(7-7')(^e-‘^)™(:e-i^’'*(“si-7')nasi-7').). 

Considering only thermal fluctuations, we can suppose 
that the LO is a displaced thermal state, i.e., Phoil') = 



exp(—|7'P/n)/(7rn). This allows us to compute the 
Gaussian integral (38) as 




1 


1 + An 


^:exp 


A 


1 + nA 


K7) :), 


where A = rjtm/N S [0,1]. Note that the prefactor of the 
initial exponent A in Eq. (37) is reduced to A/(l + hA) in 
Eq. (39). This can be regarded as a diminished quantum 


efficiency. The analytical formula (39) may be used to 


estimate the impact of thermal LO fluctuations on the 
nonclassicality probes. 

An example of the impact of a thermal LO on the mea¬ 
surement of a nonclassicality probe can be seen in Fig. 
For the even coherent state the nonclassicality cannot be 
determined for strong fluctuations of the LO. Namely, in 
the studied scenario the LO is much stronger than the SI, 
|ap ^ |/3p, which means that the click detector mainly 
detects the noisy LO. Thus, experiments should avoid 
such intensity relations, or the LO fluctuations have to 
be minimized. 





0.50 


0.00 


FIG. 8 : (Color online) The nonclassicality criterion (22 1 is 
shown for the even coherent state |a+) snperimposed with a 
LO, where a = 1 and /3 = 4. The negative normally ordered 
variance of 7 i((p) is analyzed depending on the thermal photon 
number h of the LO and the phase </ 9 . The detection and beam 
splitter parameters are 7 = 1, iz = 0, A = 8 , t = 4/5, and 
r = 3/5. 


A second effect, which is due to the LO, is a mode mis¬ 
match. This imperfection occurs whenever the modes of 
the LO and the SI field do not have a perfect overlap 
when combining them on a beam splitter (see Fig. [^. 
The rigorous derivation of the resulting modifications 
based on the spectral response of the detector can be 
found in Appendix [B| 

Here, let us discuss solely the results. The moments 
change to 


introducing an additional noise term v. In particular, 
this noise contribution is proportional to the intensity of 
the LO beam. 


V oc | 7 |^ 


(41) 


(see Appendix E- As we have seen before [Eq. ([M|)], 
such a noise contribution will not affect the nonclassical 
properties of the matrix of click-counting moments. This 
situation changes in combination with LO fluctuations, 
since it will contribute to the convolution (38). Therein, 
v has to be replaced hy v + v. 


V. MULTI-PORT HOMODYNE DETECTION 



FIG. 9: (Color online) Balanced eight-port homodyning 

scheme consisting of four 50:50 beam splitter and a tt/2 phase 
shifter. The SI and LO are each fed into an input port, 
whereas two input ports, vac.l and vac. 2 , remain unused. 
The outgoing fields are detected with click-counting detectors. 
Nonlinear click quadratures can be subsequently obtained by 
subtracting the click counts of the detectors in positions 1 and 
2 as well as positions 3 and 4. 


After this detailed consideration of imperfections, let 
us proceed by generalizing the formalism of four-port de¬ 
tection with click-counting detectors to prominent bal¬ 
anced multi-port homodyning schemes. A standard bal¬ 
anced setup is the well-known eight-port scheme, shown 
in Fig.j^for four click detectors. Its input-output relation 
can be written as [5] 


/ai\ 


(-1 1 

1 

02 

1 

I 1 

-1 

03 1 

~7t 

I i 

1 

XdJ 

\ 1 -* 

1 


/ asi 

olo 

^vac. 1 

V^vac.2 


(42) 


{:[N - :^(v?)]'":) mismatch 


( 40 ) 


Again, we assume that all click detectors have the same 
characteristics. Thus, generalizing the previously intro- 
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duced approaches, we get 
(:(di - d-2)''(d3 - d-4)':) 

X sinh* 


(43) 


v\P\ 

AN 


x{ip) 


»7|/3U/ , tt, 

4iv"^<‘^+2> 


:>. 


with k,l = 0,..., N. The conjugate momentum to x{ip) 
is p{ip) = x{(p + |). This can be adopted to define the 
nonlinear momentum operator, 


P{(p) =2iVe"l3?^l'^l"+‘^':e"a?f”sinh 


V\P\ 


(44) 


for the nonlinear quadrature operator, here in the form 


X{<f) =2Ne ■'■‘^1:6 ^w^-sinh 


V\P\ ^ 


2 (45) 


Let us point out that X{(p) in Eq. (28 1 for the four-port 


scheme includes terms which scale with rj/{2N), whereas 
for the eight-port scheme we have a scaling with 'q/{AN) 
[see Eq. (45)]. This deficiency could be corrected by tak¬ 
ing half the numbers of APDs for the eight-port homo¬ 
dyne detection. Using photoelectric detectors, such a 
scaling also occurs 0163]. There, however, the correc¬ 
tion requires us to have a doubled quantum efficiency in 
the eight-port scheme instead, which is a quite demand¬ 
ing task. In click detection, the total efficiency of the 
detector system, consisting of N APDs, each having a 
quantum efficiency 77 , is given by the fraction p/N. This 
allows us to modify N 01 rj for manipulating the overall 
efficiency. This relates to the findings in Ref. where 
it has been demonstrated that it can be advantageous to 
have fewer on-off detectors in some scenarios. 

Finally, we can formulate nonclassicality in terms of 
variances, 0 > (:[AA((^)]^:) or 0 > (:[AP((/ 3 )]^:), or we 
can uncover nonclassical correlations between click posi¬ 
tion and momentum via 

0 > {■\xx{p)f-){:\xP{p)Y-) - {■.xx{p)xP{p)-)\ 

(46) 

In this form, the nonclassicality condition relates to a vi¬ 
olation of a normally ordered version of the Schrddinger- 
Robertson uncertainty relation j^. Let us stress again 
that not only can moment-based nonclassicality crite¬ 
ria be constructed from the second-order difference mo¬ 
ments, but more general criteria may be considered as 
well. For instance, the variance of the sum of click events 
from different detectors could be used to certify the quan¬ 
tum character of the odd coherent state. 

Examples for the single quadrature variances and the 
covariance are given in Fig. Nonclassicality is de¬ 
termined for the even coherent state. Here, it is worth 



(:[AX(^)]":>(:[AP(^)]2:>-<:AX(^)AP(^):>" 


-9.6560 

-9.6562 

-9.6564 


0 


27r 


FIG. 10: (Color online) Both click quadratures variances [top 
left: X{ip), top right: P(v3)] and the covariance in Eq. | |46| ) 
(bottom, scaled xlO^) are shown as a function of the phase 
ip for the even coherent state |q:+) with a = 1. As considered 
before, the detection parameters are r] = 50%, p = 0, A = 8, 
and the LO amplitude is |/?| = 4. 


p. The nonlinear position and momentum variances are 
negative only for small phase intervals. 

In general, the method derived here applies to all kinds 
of multi-port homodyning schemes. For example, a six- 
port scheme with photoelectric detectors was proposed 
in Ref. [63| and could be similarly formulated with click¬ 
counting techniques. Since the six-port scenario includes 
only one vacuum input, it yields smaller attenuations 
than the eight-port scheme with two vacuum inputs dis¬ 
cussed here (see Fig.[^. 


VI. MULTI-MODE MEASUREMENTS 

So far, we have considered only a single SI field. Here, 
we outline the multimode scenario of phase-sensitive click 
detectors at different positions, which is especially in¬ 
teresting for determining quantum correlations between 
beams. For simplicity, we restrict our discussion to two¬ 
mode balanced homodyne detectors as one example for a 
correlation measurement (see Fig. 11). A generalization 


to multi-port schemes with more than two spatial modes 
or with other homodyning schemes is straightforward. 

In Fig. 0 each of the signal modes asi.i and asi .2 is 
fed into a balanced homodyne. Therein they are super¬ 
imposed with the corresponding LO modes of coherent 
states |/3i), with i = 1,2. The click difference counts in 
each balanced detectors yiel d two nonlinear quadratures 


pointing out that the evaluated condition (46) includes Xi{pi) and X 2 ( 722 ) (see Sec. HIB). Single-mode nonclas- 


very small oscillations, and it is negative for any phase sicality for the SI in mode i = 1 or f = 2 is verified if the 
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FIG. 11: (Color online) A two-mode correlation measure¬ 
ment scheme consisting of four click counters. Each of the two 
SI and LO beams are fed into a balanced homodyne detection 
(“HDl” and “HD2”) four-port scheme. The individual, non¬ 
linear click quadratures of the two modes, Xi{ipi) for i = 1, 2, 
are subsequently obtained by subtracting the click counts in 
each setting. 




'Pi Pi 





corresponding normally ordered variance is negative, 

0 > (: :)• (47) 

Additionally, nonclassical two-mode correlation between 
the two signals can be inferred from a negative covari¬ 
ance: 


0>(:[A1i(:^i)]^:)(:[AX2(:^2)]2:) 

-(:AAi(^i)Al2(:/92):)", 


(48) 


which yields a nonlinear two-mode squeezing. It is also 
worth mentioning that a quantum correlation between 
the Sis in terms of a sum of the local balanced homodyne 
click counting devices can be formulated, as done for the 


single-mode case in Sec. IIIB 


Let us apply this approach. The natural extensions of 
the single-mode even and odd coherent states in Eq. i> 
are the two-mode odd and even coherent states. 


FIG. 12: (Color online) The click quadrature variances (471 
(top left: i = 1, top right: i = 2) and the covariance in 
Eq. (481 (bottom, scaled xlO^) are shown as a function of 
the phases (pi, i = 1, 2. Nonclassicality of the two-mode even 
coherent state |a+), with a = 1, is revealed for click-counting 
detection parameters rj = 50%, = 0, A = 8, and the LO 

amplitudes are |/3i| = |/? 2 | =4. 


that each APD produces not more than one click in a 
given measurement interval. The corresponding nonclas¬ 
sicality conditions, e.g., for the click counterpart of the 
photon antibunching effect have also been for¬ 

mulated in terms of normally and time-ordered matrices 
of click-counting moments. 


VII. SUMMARY AND CONCLUSIONS 


^(2)^ ^ Iq, q;) ± I - a,-a) 
^ v/2[l±exp(-4|ap)] 


(49) 


The results for the second-order criteria are illustrated 
in Fig. |12[ For the given parameter range, the single¬ 
mode variances of the two-mode even coherent state are 
negative for a small neighborhood of 7r/2 and 37r/2. In 
contrast, the nonclassical covariance between the modes 
is negative for almost all phase values even for a detection 
efficiency of 50%. 

Let us remark that multi-time, multi-detector cor¬ 
relation measurements as reported in Ref. [64] can be 
achieved with a similar approach. The multi-time click¬ 
counting theory was established [52] under the constraint 


In summary, we have formulated a theoretical model 
for the implementation of click detectors in phase- 
sensitive homodyne measurement schemes. Based on 
such setups and employing the click-counting theory, we 
studied the verification of nonclassical light for realis¬ 
tic detection processes. Since our results have been ex¬ 
pressed in terms of analytical formulas, they may be 
helpful for predicting the experimental results and for 
formulating bounds to imperfections in realized setups. 
In addition, the moments for the considered nonclassi¬ 
cality probes can be directly derived from measured click 
statistics, making our findings easily accessible for exper¬ 
imental implementations. 

We formulated the four-port homodyning including the 

















11 


balanced and unbalanced detection scenarios. Our ap¬ 
proach in terms of click-counting detectors was compared 
with the traditional scenario using photoelectric detec¬ 
tion models. Perturbations stemming from the imperfect 
click-counting detectors, e.g., efficiencies and saturation 
effects, or the impurities of the local oscillator have been 
studied. Moreover, the influence of a mode mismatch be¬ 
tween the signal and the local oscillator has been shown 
to result in an additional dark count rate which is propor¬ 
tional to the intensity of the local oscillator. In the case 
of multi-port homodyning, we identified a nonclassicality 
condition in the form of an uncertainty relation between 
the nonlinear position and momentum operators. In the 
same fashion, quantum correlations between multiple sig¬ 
nal fields have been studied. It is worth mentioning that 
we focused our consideration of nonclassicality probes on 
second-order criteria. The extension to higher-order mo¬ 
ments was also discussed. 

We can conclude that click detectors are capable of de¬ 
termining nonclassicality in various phase-sensitive opti¬ 
cal measurement scenarios. This supports the assump¬ 
tion that click-counting devices can be employed when¬ 
ever photoelectric detectors are not available. Hence, our 
technique offers a useful set of tools for current and fu¬ 
ture experiments. It may also be the starting point for 
the development of a tomographic state reconstruction 
approach with click-counting detectors, which requires 
further studies. 
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To do so, let us expand an arbitrary matrix element 
(:(a:i + y^)"':) (A4) 


=EE 


k' 


ykyk\-,%k+k' 


k=0 fc'=0 

Now, we can define a matrix T{x,y) = {tm,k)m,k with 

0 for k > m 


im,k — 


(T)^™ for k <m. 


(A5) 


The matrix T{x,y) is an upper triangular matrix with 
non zero diagonal elements, tm,m = y™ 7 ^ 0. There¬ 
fore, 7Ta;,y) is invertible. In addition, the expansion in 
Eq. (A41 proves the transformation 


M{x,y) = Tix,y)MiO,l)T\x,y). 


(A 6 ) 


Due to this transformat ion property and because 
T{x,y)~^ exists, the claim (A3) holds true. 


Appendix B: Mode Mismatch 

We will derive the description of a mode mismatch 
between LO and SI for click detectors. This approach 
is based on multimode detection with spectral response 
functions (see [H ). A single APD is properly char¬ 
acterized by two operators for the click, Pon = 1 — 
:exp[—F]:, and no-click event, PqS = :exp[—F]:, where 
F is the detector response. Restricting ourselves to spec¬ 
tral properties (spatial and polarization degrees of free¬ 
dom can be treated similarly) we have 


Appendix A: Matrix of Moments Expansion 

Nonclassicality conditions can be written in terms of 
matrices of moments. The typically considered matrix of 
normally ordered click-counting moments reads 

M= ((:7r’"+'"':)) . (Al) 

\ / m,m' 

Note that (: 7 r:)> 0 . Let us formulate two equivalent rep¬ 
resentations of matrices: (i) in terms of powers of A — tt, 
A4’^=((:(A— 7 r)"*+™ terms of centered 

moments Att, Af''=((:(A 7 r)'"+'" these no¬ 

tations may be rewritten in the general form 

M(x,y) = ((:(xi + y^)™+™':)) , (A2) 

with x, y € M and y 7 ^ 0. Thus, we have M'=M{N, —1) 
and Af"=A4’(—(: 7 r:), 1). We will prove that the non¬ 
negativity is preserved for all real parameters x and 
y 7 ^ 0 ; that is, we claim that 

M = M(0,1) > 0 M(x,y) > 0. (A3) 


f = y dw G{uj)a^{uS)a{u}), (Bl) 

with G{uj) > 0 being the so-called spectral response func¬ 
tion of a broad band detector [2]. The annihilation op¬ 
erator for a frequency w at the detector is d{uj). For 
simplicity, we assume a vanishing dark count rate for the 
time being. 

It has been derived in Ref. [S2] that the click statis¬ 
tics for multimode fields may be expanded in terms of 
expectations values 

/(A) = (:exp[-Af]:), with A G {0/A,..., A/A}. (B2) 

Moreover, the semi-classical representation in terms of 
the Glauber-Sudarshan representation of the SI field, 
Psi = / d?Oi P{a)\a) {a\, allows one to consider coherent 
SI first, la) = |a)si, and generalize the considerations to 
arbitrary states by integrating the result with the P func¬ 
tion. The spectral decomposition of the LO and SI is 

Oi = y du!fi{oj)di{uj), foriG{LO,SI} (B3) 
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and f duj\fi{oj)\‘^ = 1. Thus, the coherent SI and LO 
states can be written as 

|a)si = |/si(w)a)c,; and |/3)lo = \fho{uj)P)u, 

UJ LO 

(B4) 

respectively. Applying a beam-splitter transformation to 
the spectral modes [1 = |t(w)p-|-|r(tLi)p], we have for one 
output port of the beam splitter 


a{uj) = t{uj)dsi{io) + r(a;)dLo(w). (B5) 


Hence, the desired expectation value (B2) reads 


Now we can define the overall quantum effi¬ 
ciency ?7t=(t/sii i/si); the coherent displacement 
7 = -[{tfsi, rfLo)/{tfsi, tfsi)]l3, and a rate for the mode 
mismatch 


(t/spt/si)(r/Lo,r/Lo) - |(t/si,r/Lo)| 
(tfshtfsi) 


W?- (B9) 


Note that Cauchy-Schwartz inequality implies v >Q. Us¬ 
ing the P representation and the normal-ordering pre¬ 
scription, we get our desired quantity (B2) for arbitrary 
SI states psi from Eq. (B81 as 


I (A) = exp [-X{tfsia + rf^oP, tfsia + c/lo^)] , (B6 ) 
using the inner product 

{a,b) = j doj G{uj)a* (B7) 


as 


The derived expression (B61 can be further rewritten 

(B8) 


/(A) = exp 
X exp 




, («/si,r/L0) o 

+ (t/si.t/si) P 


jy (t/si7/si)(»~/LO,'r/Lo)-|(t/si,'r/Lo)| I o|2 
(tfsiytfsi) P' 


/(A) = (: exp [-A?7t(osi - 7)'''(asi - 7 ) - ■)■ (BIO) 


In comparison to a perfect mode matching (t/si = c/lo 
=> iz = 0), we have (i) the same structure of an expo¬ 
nential of a displaced photon number operator h( 7 ) [see 
Eq. @] and (ii) an additional dark count rate v which 
is proportional to the intensity of the LO [see Eq. ( |B9[ )] , 
or, equivalently. 


i> oc jyp. 


(Bll) 
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